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Abstract 



Multidimensional cosmological model describing the evolution of a fluid with shear and bulk 
viscosity in n Ricci-flat spaces is investigated. The barotropic equation of state for the 
density and the pressure in each space is assumed. The second equation of state is chosen 
in the form when the bulk and the shear viscosity coefficients are inversely proportional to 
the volume of the Universe. The integrability of Einstein equations reads as a colinearity 
constraint between vectors which are related to constant parameters in the first and second 
equations of state. We give exact solutions in a Kasner-like form. The processes of dynamical 
compactification and the entropy production are discussed. The non-singular D-dimensional 
isotropic viscous solution is singled out. 



g^; PACS numbers: 04.20.J, 04.60.+n, 03.65. Ge 

> : 

X 1 Introduction 

a' 

Last two decades have witnessed an increase of interest to the multidimensional cosmol- 
ogy (see, for instance, [5-7,9,11-18,20-25,27-29,31,32,39,41,42,45-48] and references therein). 
According to this theory one assumes that the Universe had a higher dimension at a very 
early stage of its evolution, and that quantum processes have been responsible for the (topo- 
logical) partition of the space, which provides us at present with the usual 3-dimensional 
(external) space, in addition to internal space(s). The manifold which accounts for such a 
multidimensional spacetime has the following topology 

M = Rx M 1 x . . . x M n , (1.1) 

where R stands for the cosmic time axis, and the product with one part of manifolds 
Mi, . . . ,M n gives the external space, when the remaining part stands for internal spaces. 
The classical stage of the evolution is governed by the multidimensional version of Ein- 
stein's equations. According to the classical description, the Riemann curvatures of spaces 
Mi, . . . , M n are assumed to be constant (Einstein spaces). Such a model is the simplest mul- 
tidimensional generalization of the space-time upon which the Friedman-Robertson- Walker 
(FRW) world model is based on. One easily understands that the use of extra dimensions 
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(for the physical space-time) can be a sensible scenario only at primordial epochs, since 
the standard FRW world model is known to be in a sufficiently good agreement with the 
observational constraints down to a quite primordial epoch such as the nucleosynthesis era. 
Hence, it is clear that a reduction process (called dynamical compactification of additional di- 
mensions) is required before such an epoch, to make the internal spaces contracting themself 
down to unobservable sizes. 

Herein, we assume that the cosmic fluid (the source of the gravitational field at early 
stages) is viscous, which might simulate high energy physics processes (such as the particles 
creation). The effects related to viscosity in 4-dimensional Universe were studied through dif- 
ferent viewpoints (see e.g., [2-4,8,10,26,30,33-38,44]). Before developing the 
multidimensional model let us briefly discuss (extensive review of the subject was given 
by Gron [19]) the main trends in 4-dimensional cosmology with viscous fluid as a source. 
First, Misner [3~3] considered neutrino viscosity as a mechanism of reducing the anisotropy 



in the Early Universe. Stewart |4(| and Collins and Stewart [T(| proved that it is possible 
only if initial anisotropies are small enough. Another series of papers which concerns the pro- 
duction of entropy in the viscous Universe was started by Weinberg | 4"3}| . Both isotropization 
and production of entropy during lepton era in models of Bianchi types I,V were considered 
by Klimek P5| . Caderni and Fabbri M calculated coefficients of shear and bulk viscosity in 
plasma and lepton eras within the model of Bianchi type I. The next approach is connected 
with obtaining singularity free viscous solutions. The first nonsingular solution was obtained 
by Murphy |}5] within the flat FRW model with fluid possessing a bulk viscosity. However, 
Belinsky and Khalatnikov [Q, [3J showed that this solution corresponds to the very peculiar 
choice of parameters and is unstable with respect to the anisotropy perturbations. Other 
nonsingular solutions with bulk viscosity were obtained by Novello and Araujo H36| , Romero 



38| , Oliveira and Salim j37|. 

The crucial feature of each viscous cosmological model is assuming of the so called "second 
equation of state", which provides us with the viscosity coefficients dependence on time. 
Further we denote by ( and rj the bulk and shear viscosity coefficients, correspondingly. 
Murphy integrated the 4-dimensional flat FRW model with bulk viscosity by assuming 
( ~ p (as second equation of state), where p is the density of the viscous fluid. Belinsky and 
Khalatnikov |3[] studied the behavior of this model as well as homogeneous anisotropic 
models of Bianchi types I and IX by means of qualitative methods with more general second 
equations of state (, r\ ~ p u , where v is constant. Lukacs [^] integrated the homogeneous 
and isotropic 4-dimensional model with a viscous pressureless fluid and a second equation 
of state given by ( ~ [scale factor] -1 . A curvature-dependent bulk viscosity was studied in 



multidimensional cosmology by Wolf f4"?fl . Recently, Motta and Tomimura [33J] studied a 
4-dimensional inhomogeneous cosmology with a bulk viscosity coefficient which depends on 
the metric. In our previous papers |T5| , |TT|| exact solutions in multidimensional models with 
bulk viscosity were obtained and their properties were studied for the following type of the 
second equation of state: ( ~ [ volume of the Universe] -1 . 

The aim of the present investigation is to integrate the Einstein equations for a 
multidimensional cosmological model formed by a chain of Ricci-flat spaces and a cosmic 
fluid possessing both shear and bulk viscosity. The second equations of state are chosen 
in the following form of metrical dependence of the bulk and shear viscosity coefficients: 
(,rj ~ [volume of the Universe] _1 . 
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The paper is organized as follows. In Sec. 2 we describe the model and get basic 
equations. To integrate them, we develop some vector formalism proposed in our previous 



papers |T4|, [16[. Thermodynamical concepts in multidimensional cosmologies are defined 
in Sec. 3, where a formula which provides us with the variation rate of the entropy is 
derived. The equations of motion for the special set of parameters in the first and the 
second equations of state are integrated in Sec. 4. Exact solutions are given in a Kasner-like 
form, their physical properties are investigated in Sec. 5, where the process of dynamical 
compactification and the entropy production are also defined. 



2 The model 

Let us have the following metric 

n 

ds 2 = -e 2 ^dt 2 + ]Texp[2x*(t)]d^, (2.1) 

i=l 

on the manifold defined in Eq. ( |1 . 1| ) , where dsf is the metric of the Einstein space Mi, ^(t) 
and x l (t) are scalar functions of the cosmic time t. The dimension of this manifold is given 
by D = 1 + Yn=i-Ni, where Ni = dimMj. Herein, for reason of simplicity, only Ricci-flat 
spaces Mi, . . . , M n are assumed (i.e., the components of the Ricci tensor for the metrics dsf 
are zero). One easily obtains the non-zero components of the Ricci-tensor for the metric 
defined in Eq. ([H]) (see |3|): 

R°o = e-^*)^Ai(i') a +7o-77ol, (2-2) 

K; = e- 2 ^(^ + ( 7o -7)iOC' (2-3) 

for % = 1, . . . , n, where the indices m 8 and rii run from D — J2]=i Nj to D — Z)" =i Nj + and 
7o = £?=iA^. 

A viscous fluid is characterized by a density p, a pressure p, a bulk viscosity coefficient (, 
a shear viscosity coefficient i], so that the (standard form of the) energy-momentum tensor 
reads 

7? = pu?u v + (p - C0)PH - 2 W », (2.4) 

where is the D-dimensional velocity of the fluid, 9 = denotes the scalar expansion, 
Py = $1/ + y^u v is the projector on the (D — l)-dimensional space orthogonal to w M , and 
a v — \ ( u a;/3 + U/3-a) P^Pf? — (D — 1) _1 ^P^ is the traceless shear tensor (defined as usual). 

By choosing the D-dimensional velocity so that = S^e" 1 ^ (the comoving observer 
condition), we obtain 



9 = 7o e-^, (2.5) 

{u»u v ) = diag(-l,0,...,0), (2.6) 

(P?) = diag(0,l,...,l), (2.7) 

K) = e-^diag^,^ 1 -^)^^...,^--^)^;), (2.8) 
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where ki, U = l,...,Ni for i = 1, ...,n. The function / y(t) determines a time gauge (a 
harmonic time gauge for y(t) = 70 and a synchronous time gauge for j(t) = 0), see Eq. Q2.1|); 
note that the harmonic time t and the synchronous time t s are related by dt s = exp[7 ]dt. 
By assuming anisotropy properties for the pressure and the bulk viscosity, with respect to 
the whole space Mi x . . . x M n , one has 

(T^)=diag(-p,rf; i ,... ) ^f;), (2-9) 

where 



Pi = Pi - e 7 



C*7o • --M.;-' 



and pi, resp. is the pressure, resp. the bulk viscosity coefficient, in the space described 
by the manifold Mj. Furthermore, we assume that the barotropic equations of state hold 

p, = (l-h t )p(t), (2.11) 

where the hi are constants. One easily shows that the form of the equation of motion 
(VmTq 1 = 0) for a viscous fluid described by a tensor given by Eq. fl2.9| ), is given by 

n 

p + J2N l x l (p + p*) = 0. (2.12) 

i=l 

The Einstein equations R% — ^S^R = k 2 T^, where n 2 is the gravitational constant, can 
be written as = k 2 (T^ — j^Stf). Further, by using the equations R® — ^SqR = k 2 Tq and 
R%i = n 2 {T^ - ^5^), Eqs-d^SED g ive the following equations of motion 

n 

^iy,(x 4 ) 2 -7o 2 = -2K 2 e>, (2.13) 

i=l 

e 7 p -hi + k ~ x ' 



+ (70 - = t^e 1 

+ 7o ( -d + 



D-2 

D-2 



2^ V [x i -^- [ ). (2.14) 



We use an integration procedure which is based on the n-dimensional real vector space 
R n . Let ei, . . . , e n be the canonical basis in R n (i.e. e\ — (1, 0, . . . , 0) etc. . . ), and (, ) denote 
a symmetrical bilinear form defined on R n by 

(d, ej ) = 8 ij N j - NiNj = G i5 . (2.15) 

It has been used as a mini-super-space metric for cosmological models (see [20-25]). Such 
a form is non-generate and has the pseudo-Euclidean signature (— , +,..., +). With this in 
mind, a vector y e R n is time-like, resp. space-like or isotropic, if (y, y) takes negative, resp. 
positive or null values; and two vectors y and z are orthogonal if (y, z) = 0. Hereafter, we 
use the following vectors 

x = x l ei + . . . + x n e n , (2.16) 
u = u 1 e 1 + ... + u n e n , tt < = fr- * » ^ = N i^ ( 2 - 17 ) 

£ = ^ ei + ... + Ce n , e = d- » & = JV &. ( 2 - 18 ) 
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where covariant coordinates of the vectors are introduced by the usual way. Moreover, let 
us denote u d the particular vector given by Eq. (2.17) with h i= \^,^ n = 1 (it is related to dust 
in the whole space, see Eq. (|2.11 )). One has 

(u d )i = Ni, u d = — j - , (u d ,u d ) = ~ _ 1 (u d ,x) = j . (2.19) 

Thus, using Eqs. (2.16-2.19) we rewrite the Einstein equations (2. 13), (2. 14) in the form 

(x,x) = -2/? 2 e 27 p, (2.20) 
x + ({u d , x) — 7 + 2r]K 2 e' 1 ^ x = 



-u 



- K 2 e?(u d ,x) U - i ^ u s j u *\ > ( 2 ' 21 ) 

where the formal dependence on p in Eq. (2.21) has been canceled, according to Eq. (2.20). 
Moreover, Eq. (gj^) can be written as 

p + (2u d - u, x)p - e" 7 Ur){x, x) + (u d , x) {£ - , x^j = 0. (2.22) 

To integrate Eq.(2.21) one needs second equations of state, involving the bulk viscosity coef- 
ficients Cij • • • > Cn an d the shear viscosity coefficient rj. Let us assume that these coefficients 
are proportional to exp[— 70] (or inversely proportional to the volume of the Universe), i.e. 

77, Ci ~ [scale factor of M!]- dim(Ml) • . . . • [scale factor of M n ]" dim(Mn) , (2.23) 

which means (from a physical viewpoint) that the viscosity decreases when the space Mi x 
. . . x M n expands. The integrability of the basic equation (provided the second equations 
of state) is ensured when the vectors u, £, u d are either colinear or orthogonal (with respect 
to the mini-super-space metric) in some combination [15,18]. Herein, we suppose that these 
vectors are colinear, which means that the viscous fluid has identical properties in the internal 
space(s) and the external space. Hence, all these assumptions, for the pressures and the 
viscosity coefficients, allow us to write 

Pi = (1 — h)p or u — hu d , (2.24) 
d = ^e~^ = ( or ^ = Cu d =^e-^u d , (i = l,...,n), (2.25) 



Ik- 

where ( , r/ and h are constants. 



V = ^\ (2-26) 



3 Multidimensional Thermodynamics 

According to [43,48], let us summarize thermodynamics principles in such a 
multidimensional cosmology. The first law of thermodynamics reads 

TdS = d(pV) + Vj2^, (3.1) 

i=l V i 
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where Vi stands for a fluid volume in the space Mj, when V — Vi ■ . . . ■ V n is a fluid volume in 
the whole space, and S is an entropy in the volume V. By assuming that the baryon particle 
number in the volume V is conserved, Eq. ( |3.1| ) transforms to 



nTs = p + pY^ + Pi N iZ 



X 



(3.2) 



i=l 



i=l 



where s = S/N-q, resp. n = Nb/V, stands for the entropy per baryon, resp. the baryon 
number density. Let us remind that exp fx*] is the scale factor of the space Mj (of dimension 
Ni). For a perfect fluid (Q = 0, r\ = 0), the comparison between Eqs. (|2.12j , |3.2|) gives the 
entropy conservation ( i.e., s is constant). Similarly, we obtain also the temperature (see 



48]). From Eq. (O) we have 



' dp ~ 



-pNi - Pi Ni = (hi - 2)N lP , (j ± i), 



and then 



p = K(s) exp 



x 



i=l 



(3.3) 



(3.4) 



where K(s) is an unknown function (which reads in term of the entropy per baryon s). Using 
Eqs. (3.2), (3.4) we get 



— - ] = nT = K'(s) exp 
os / _ t 



E( 2 - W 



,i=i 



(3.5) 



For a perfect fluid, we have K'(s) = 1/B where B is a constant, then 



nT = — exp 



J2(hi - 2)N l 



x 



i=l 



— exp[(w - 2u d ,x)}. 

JD 



(3.6) 



For a fluid with a bulk and shear viscosity, the comparison between Eqs. ( |2.12| , |3~2D pro- 
vides us with 

nTs = e" 7 ^2r){±, x) + (u d , x) (£ - 2 — ^-^ u d , x)j . (3.7) 

Such a formula gives the variation rate of entropy per baryon in multidimensional cosmology 
on the manifold M = R x Mi x . . . x M n with viscosity. The entropy production can be 
calculated if the temperature of the fluid is known. Herein, we suppose that the temperature 
is given by the perfect fluid formula Eq. (3.6), which is valid with sufficient accuracy when 
effects of viscosity are small. Hence, Eqs. (3. 6), (3. 7) give 



s = Be {2ud - U ' x) ^ 2 V {x, x) + (u d , x){£-2 



2wd 

(Ud, Ud) 



,x) 



(3.8) 
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4 Exact solutions 



According to assumptions given in Eqs. (2.24-2.26), the basic vector equation, see Eq. (2.21] 
reads 



x + {{udi x) — 7 + 2r/K 2 e 7 ) x 



^(x, x) - K 2 e'(u d , x) ( C - . 2V . 
2 \ {u d ,u d ) 



In order to integrate such a (vector) equation, we use the orthogonal basis 

7 r> h, ■ ■ ■ , fn £ -R n ; 

[Ud, Ud) 

where the orthogonality property reads 

fad, fj) = 0, (fj, fk) = S jk , (j, k = 2, . . . , n). 



(4.1) 



(4.2) 



(4.3) 



Let us note that the basis vectors / 2 , . . . , f n are space-like, since they are orthogonal to the 
time-like vector u d . The vector x G R n decomposes as follows 



x = (u d ,x) 



u d 



(4.4) 



Hence, the basic vector equation given in Eq. (|4.1|) reads in term of coordinates in such a 
basis as follows 

(fj, x) + ((u d , x) - 7 + 2r^V) (fj, x) = (4.5) 



(u d ,x) ( (u d ,x) -7 2 \ /l 



+ 



J'=2 



(4.6) 



which provides us with a set of equations (for j = 2, . . . ,n). For the metric dependence of 
the viscosity coefficients, see Eqs. (2. 25), (2. 26) Eqs. (|4.5|,|4.6|) read 



(fj,x) + Vo(fj, x) = 0, (j = 2, . . . , n) 



\Ud,x} - ^(u d ,x) 2 + (u d ,u d ) Uo(ud,x) - ^Yl(fj,x) 2 ] = 



in the harmonic time gauge 



7 = 7o = (ud, x) = ^Ni 



x 



(4.7) 
(4.8) 

(4.9) 



i=l 



Such a set of equations is integrable for any values of the constant parameters h, r/o and Co- 
Let us first assume models with 

h^O. (4.10) 



The integration of Eq. (|4.7f ) gives 

(fj,x) 



fpj _)_ gi if r/o = 0, 
e-*V + g-? if r/o ^ 0, 



(4.11) 



where p> and g J are arbitrary constants. The Kasner-like form solution can be written in 
term of vectors a, f3 G R n , defined as follows 



a = E p'fi = E a%e ^ P = E tffi = E 

j=2 i=l j=2 i=l 



(4.12) 



where a* and /3 l are their coordinates in the canonical basis ei, . . . , e n . By using the orthog- 
onality conditions, we obtain 



(a, u d ) = £ aW, = 0, (/?, u d ) = £ ^ = 0, 

i=l i=l 
n 2 n 2 



8=1 



i=2 



3=2 



(4.13) 

(4.14) 
(4.15) 



t=i 



where the constants a* and /3* may be called Kasner-like parameters, because of the existence 
of these constraints. By using Eqs. ( f4.4| , |4.11| . f4.12| ), we obtain 



where the function 



x = (u d ,x) 



a(t) 



Ud 



(Ud, Ud) 



+ a(t)a + (3, 



t if r]o = 0, 
e"** if 77o ^ 0, 



(4.16) 



(4.17) 



By substituting the functions (fj,x) into Eq. (O) we obtain the following equation for the 
unknown function (ud, x) 



h h 
(u d ,x) - -(u d ,x) 2 + (u d ,u d )( (u d ,x) = --A 2 a 2 (t), 



where 



.4 



IP-1 
D-2 



(a, a). 



(4.18) 



(4.19) 



The equation Eq. (|4.18| ) has been integrated for the non- viscous model r] = ( = (see, e.g., 
14J ) and for the model with Co 7^ and 770 = in [15, 17]. For the model with rjo 7^ 0, it can 



be reduced to the modified Bessel equation 



,d 2 z dz 



+ T- ( T ' + V Z )Z = 



dr 2 dr 



by means of the transformation 



r = -\h\Ae- vot , 
2 \ 1 

(u d ,x) = ^r— h\\T~ v z{r) 
h dr 



(4.20) 

(4.21) 
(4.22) 
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in the non-trivial case (a, a) 7^ 0, where the constant 

_ D - 1 Co 
27 D — 2 2i] ' 

The general solution of the modified Bessel equation is given by 



(4.23) 



(4.24) 



see where I\ v \{r), resp. K\ v \(t), is the related modified Bessel function, resp. Mac-Donald 
function. 

Finally, the results of Eqs. ( [4.7| , |4.8|) integration for various constants Co an d f]o can be 
presented in Tab. 1 

Table 1 





Co = 


Co 7^0 


770 = 


Solution II with (a, a) 7^ and 
Solution I 


Solution II 




Solution III with (a, a) 7^ and 
Solution I 


Solution III with (a, a) 7^ and 
Solution II with (a, a) = 



The solutions I, II, III in term of scale factors are the following 

• Solution I : 

^ = ^\C l + C 2 t\- 2/[h{D - 1)] . 

• Solution II : 



Q Q -(A-Co)ht/2 + q (A+Co)ht/2 



-2/[h(D-l)] 



where 



<b=£^. A 



D-2 h 



D-l (D-lCv 



Solution III : 



exp (a l e-^ + fi) (t~ v \C x I\ v \{t) + C 2 K H (r 



-2/[h{D-X)} 



(4.25) 



(4.26) 



A 2 + (1 (4.27) 



(4.28) 



where the variable r > is given in Eq. ([4.21 ), the constants A in Eq. (|4.19|) and v in 
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In formulas (4. 25), (4. 27), (4. 28) Ci-x$ are integration constants such that Cf + Cf > and 
i — 1, . . . , n. The Kasner-like parameters a 1 and f3 l obey the relations given in Eq. (4.13). 
The set of equations given in Eqs. (|4.7|,f4.8|) is easily integrable in the case 



0. 



(4.29) 



which, with the barotropic equation of state in mind, relates to Zeldovich or stiff matter. 
The results are given as follows : 



Solution IV : for i 



n 



exp 



a*a(t) + ft 



exp (Ci + C 2 t) if Co = 

exp(C 1 + C 2 exp(Co§^)) if Co t^O, (430) 



where 6*1=1,2 are arbitrary constants, and the function a(t) is given in Eq. (|4.17|) . 



5 Discussion 

Let us remind the multidimensional generalization of the well-known Kasner solution [ 21 ], 
it reads (for the synchronous time t s ) as follows 

n 

ds 2 = -dt 2 s +J2At 2 / ds 2 . (5.1) 
i=i 

Such a metric describes the evolution of a vacuum model defined on the manifold R x Mi x 
. . . x M n , where the Mj are Ricci-flat spaces of dimension iVj with the metric ds^, A4 are 
arbitrary constants and are the Kasner parameters, which satisfy the relations 



1 1 



EJVie^l, E^(^) 2 = l- (5-2) 

i=l i=l 

The e l and the Kasner-like parameters a 1 (used in the above formulas for the exact solutions) 
are related as 

a 1 1 

e J = ±- + -^- T , (a,a)^0. (5.3) 

By using Eq. ( |4.16|) (i.e., a general decomposition of the vector x G i? n ) and the result of 
Sec. 4, see Eq. (3.8), we obtain the variation rate of entropy 

s = - 2 e~ h ^ (Co7o + a)a 2 {t)) , (5.4) 

which shows that the entropy increases when ( dt > and r)odt > 0. Further we assume 

Co>0, 7? >0, (5.5) 
so harmonic time t (as well as synchronous time t s ) increases during the evolution. 
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One easily shows that the weak energy condition T^v v v ili > 0, for any D- dimensional 
non space-like vector v v (and thus as well as for the 4-dimensional case), applied to the 
stress-energy tensor given in Eq. ( |2.9| ), can be written as inequalities 

P>0, P + p*>0, (i=l,...,n), (5.6) 

where p, resp. p*, is the density, resp. the effective pressure, of the fluid. 

The dominant energy condition T^v v v il > and T^T^v u v x < 0, for any non space- 
like vector v^, applied to the stress-energy tensor given in Eq. (|2.9|) , reduces to inequalities 
defined in Eq. ( |5.6|) , with the following additional condition 

p-p*>0, (i=l,...,n) (5.7) 

Notice that due to the weak energy condition the following restriction on the constant h 
(taken from the barotropic equation of state (2.11)) arises in the nonviscous case: h < 2. 
The dominant energy condition for the nonviscous stress-energy tensor implies: < h < 2. 

It is important to note that Solution I and Solutions II, IV for zero Kasner-like parameters 
(a 1 = 0, i — 1, . . . , n) are isotropic, since the spaces Mi, . . . , M n have identical scale factors. 

Further we discuss the solutions obtained, which are of interest within multidimensional 
or 4-dimensional cosmology. 



5.1 Non Viscous Models 

For a better understanding of the viscosity effect on the dynamics, we first outline the 
properties of non viscous models. 



5.1.1 



The isotropic non viscous model is described by Solution I (h 0) and Solution IV {h = 0) 
with a 1 = (for i = 1, . . . , n), which represents the multidimensional generalization of the 
flat FRW model. It is the steady-state model 



exp 



D 



-t. 



p = const 



for h = 2 and shows a power-law behavior 



(5.8) 



(5.9) 



for h ^ 2, where t s is the synchronous time (stationary solution with zero density is also 
possible); let us call it as Friedman-like behavior. 



5.1.2 

The anisotropic non viscous model for h e (0,2] is described by Solution II with C = 
and (a, a) ^ 0. One easily shows that the integration constants Ci = i 2 have to satisfy the 
condition CiC 2 < 0, otherwise p < 0, Vt. By using this condition, we obtain from Eq. (4.26) 
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-.70 



\smh[Ah(t-t )/2}\ 



-2/h 



(5.10) 



where to = can be chosen (with no loss of generality). Then, for a suitable choice of 
the integration constant of equation dt s = exp(7 )dt, one has the following correspondences 
t G (— oo, 0) t s G (0, +oo) and t G (0, +oo) <^> t s G (— oo, 0). Hence, we solely investigate 
the solution t s G (—oo, 0), since the evolution of the non- viscous fluid is reversible. 
From Eq. Q5.10 ), we obtain in the main order exp 

[7o] ~ t~ 2/h when t -> +0 (t s -> +oo), 
then one has t s ~ f 1 ~ 2 / h f or ^ g (o, 2) and t s ~ hit for h = 2. By using these relations and 
Eq. (4.26), we easily see that the multidimensional Universe shows an isotropical Friedman- 
like contraction, as defined by Eq. ( |5.9| ), in the (infinite) past. Such a conclusion is also valid 
for Zeldovich matter (h = 0). 

Let us now investigate the behavior of the non- viscous anisotropic model at t s = 0. For 
Solution II (h ^ 0), by using Eq. (5.10), we obtain in the main order exp [70] ~ exp [—At] 
when t — > +00 (t s — ► +0), then £ s ~ exp [—At]. By substituting the latter relation into 
Eqs. (4.26), we obtain in the main order 



It- 



-a>/A+l/(D-l) 



(t s -> +0). 



(5.11) 



According to Eq. (|5.3| ), the model for h G (0,2] has a Kasner-like behavior near the singu- 
larity (at t s = 0). Such a behavior describes the contraction of some spaces M\, . . . , M n and 
the expansion for the other ones. According to Eq. ( |5.2| ), the number of either contracting 
or expanding spaces depends on n (the total number of spaces) and N i= ^ n (their dimen- 
sions), but there is at least a contracting manifold and expanding one. Such dynamics is a 
mechanism of extra dimensions compactification (within the multidimensional cosmology). 

One can easily show that the non- viscous anisotropic model for Zeldovich matter (h = 0) 
described by Solution IV has almost the same properties but with the second constraint given 
Eq. fl5.2|) for the Kasner parameters substituted by Y%=i Ni(e % ) 2 = e, where e is a constant 
such that 1/(D - 1) < e < 1. 



5.2 Models with bulk viscosity 

Let us now investigate the models with bulk viscosity. 



5.2.1 

The viscous isotropic model shows interesting features. The shear viscosity is not significant 
in this case and the model is described for h 7^ by Solution II with a 1 = for i = 1, . . . ,n. 
If C1C2 < and h > then the solution can be written as follows : for i — 1, . . . , n 



Ri 1 



Pi 



D 



exp 

h exp 
1 2C 2 



D 



D-2 
D - 1 



Co(< - 1 ) 



D-2 

2N. 



Co(t ~ t ) 



-2/MD-l)] 



P, 



(5.12) 
(5.13) 



D-2 K 2 h 2 J 



exp 



i=l 



D - 1, , 

2 Co(* 

D - 2 s v 
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x I 1 — exp 
D - 1 2B( 2 



D 



D — 2 n 2 h 2 



D-2 

n 



Co(t - 1 ) 



2{2-h)/h\ 



exp 



D - 2 S0V 



-oo 



(5.14) 
(5.15) 



where we use the set of independent constants t ,Ri, . . . , R n , defined such that 

D - 1 



Ci 

c 2 



— exp 



D-2 



Co^o 



| Ci |-2/[^-l)] exp [^ 



Ri. 



(5.16) 



Let us consider this solution on the interval (— oo,£ ) for h G (0,2). Then the synchronous 
time t s changes during the evolution on the interval (— oo, +oo). Such a solution is non- 
singular and describes a monotonic isotropic expansion of the D-dimensional Universe with 
Friedman-like stage as t s — > +oo. The density p increases from zero in the infinite past to 
some maximum value and then decreases to zero at the Friedman-like stage. The maximum 
of the density is reached at t = t m = (D — 2) \n[h/2]/[( (D — 1)] + to, when p* = —p, see 
Eq. ( |5.13| ). We have p* + p < for t G (— oo, t m ) and p* + p > for t G (t m , t ), so the weak 
energy condition, see Eq. (|5.6j ), is not satisfied on the time interval (— oo,t m ). The entropy 
per baryon monotonically increases during the evolution and tends to some constant value 
in the infinite future. 

The nonsingular solution obtained by Murphy [35] within flat FRW model with bulk 
viscosity for another second equation of state: ( = constp exhibits the similar properties 
except the violating of the weak energy condition. The scale factor monotonically increases 
from zero in the infinite past and tends to the infinity at the Friedman stage of the evolution. 
The density monotonically decreases from some constant value in the infinite past and tends 
to zero according to Eq. (5.9) in the infinite future. The weak energy condition is valid for 
t s G (— oo, +oo) (but the strong energy condition p + 3p* > is not satisfied on the interval 
(— oo, t*), where t* is some constant), so from this point of view the Murphy solution is more 
attractive. 



5.2.2 



Now, let us study the properties of the anisotropic viscous model by taking into account 
only the bulk viscosity. Such model is described by Solution II with non-zero Kasner-like 
parameters a 1 and has been previously studied |L5], [L7|]- If C\Ci > in Solution II with 
a 1 7^ then the density has negative values at some stage of the evolution, which means 
that the weak energy condition is not satisfied. Hence, hereafter such solutions are not 
considered. If C1C2 < 0, then Solution II can be written as follows : for i — 1, . . . , n 



Pi 
P 



R, 



1 



svah[Ah{t - 1 )/2] 

h cosh 2 a 



-2/[h(D-i)] 



exp 



a 



1 + sinh a sinh[Ah(t — to) + a] 



(a, a) 
2k 2 



i=l 



H-Ri smh[Ah(t-t )/2] 



2(2-h)/h 



h(D 



exp[2<Cot] 




(5.17) 
(5.18) 
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(l + sinhasinhL4/i(t - t ) + a]) , (5.19) 



— — U R i ex P 



4 =1 



(O 1-7 t 



sinh 2 L4/i(t - 1 )] 



x smh 2 [Ah(t - t ) + a], (5.20) 
where we use independent constants to, Ri, ■ ■ ■ , R n defined by 

Ci = -- exp[Aht /2), C 2 = ^ exp[-i/it /2], (5.21) 

i2i = icrj-a/lM^-i)] expfyS*], (5.22) 

and a is defined by 

sinha = cosh a = — . (5.23) 

Let us note that for t G (— oo, to) the density, given by Eq. (5.19), has negative values. In 
the following, this solution is only used on the interval (to, +oo). 

If h G (0, 2) then the harmonic time interval t : (to, +oo) corresponds to the following 
synchronous time interval t s : (— oo,t°), where we choose the integration constant t° = 0. 
We can easily prove that such a model has the stage of isotropic contraction by the Friedman- 
like law defined in Eq. (|5.9| ) in the infinite past. Near the final point of the evolution t s = 
we obtain in the main order 

The final point of the evolution t s = is singular (p — > +oo as t s — > —0), and the character- 
istic of the singularity depends on the parameter (which determines the ratio between 
the viscosity parameter (o and the anisotropy parameter A). If the ratio Co /A ^> 1 then we 
obtain from Eq. (5.24) exp[x*] ~ t]/( D_1 ) near t s = 0. In such a case, the model describes the 
contraction of all spaces M 1; . . . , M n in the vicinity of the singularity. If Co/-A <C 1 then we 
obtain ^from Eq. (5.24) exp[x l ] ~ tf (i.e. the singularity is of the Kasner type). According 
to Eq. (5.20), in both cases the model describes the unbounded production of the entropy at 
the final stage of the evolution. 

Therefore, anisotropic Solution II for h G (0, 2) and CiC 2 < describes the model 
with Friedman-like isotropic contraction in the infinite past and the anisotropic Kasner-like 
behavior near the final singularity if the parameter is small enough. Under such a con- 
dition, the behavior of scale factors remains (qualitatively) the same as for the anisotropic 
non- viscous model for t s < 0. Let us note that this model satisfies the dominant energy con- 
dition during the evolution. According to Eq. (5.18), the ratio p*/p increases monotonically 
from the value (1 — h) at the Friedman-like stage in the infinite past, and tends to 1 in the 
vicinity of the final singularity. One may also consider Solution II for h > and C\ — 0, 
C*2 7^ 0. This partial solution has the similar dynamical behavior at the final stage of the 
evolution and satisfies the dominant energy condition as p* = p, Vt s . 

If the dominant energy condition is ignored then we may also consider anisotropic Solu- 
tion II for h < and CiC 2 < 0. In such a case the harmonic time interval t : (t , +oo) 
corresponds to the synchronous time interval t s : (0, +oo) . We easily show that such 
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a model has the Friedman-like singularity, defined in Eq. Q5.9Q at t s — 0. The anisotropic 
behavior is possible far from the singularity, if the parameter | £ \/A is small enough. This 
solution satisfies the weak energy condition as the ratio p*/p monotonically decreases during 
the evolution from the value (1 — h) at the Friedman-like stage, and tends to 1 in the infinite 
future. Let us note that the asymptotical behavior of this solution far from the singularity 
(t s — > +oo) is given by the Solution II for h < and C 2 = 0, C\ ^ 0. To investigate the 
possible anisotropic behavior far from the singularity let us consider this partial solution. It 
may be written in the synchronous time as following : for % — 1 , . . . , n 



e x 



i2jt «V(A+ICo|)+i/(D-i) j ts>0 , (5.25) 
P* = P = ° " t K ~°' - t~ s \ (5.26) 

. _ p-^icoi(i + ic^ +s(0); (5 _ 27) 



(D - 1)k 2 

where the constants R\, . . . , R n are such that 



n 

l[R^ = A+\Co\. (5.28) 



It is evident from Eq. (5.25) that a 1 /{A + |£o|) + 1/(D — 1) — > e l (Kasner parameter) as 
IColA^ — > 0. Therefore, if the parameter |Co|/^ is small enough, then the model describes the 
contraction of a part of the spaces Mi, . . . , M n and the expansion of another part. According 
to Eq. (5.27), the solution describes the unbounded entropy production as h < 0. 

We note, that the anisotropic viscous model described by Solution IV (h = 0) does not 
satisfy the weak energy condition as the density has negative values at some stage of the 
evolution. 



5.3 Models with bulk and shear viscosity 

Now let us investigate the model with both bulk and shear viscosity. Such a model for h ^ 
is described by Solution III with non-zero Kasner like parameters a 1 and v > 0. By using 
the following properties of the modified Bessel functions || 



I v +i(t) < I„(t), K v+ i(t) > K v {t) Vre (0, +oo), z/ > 0, 

it can be proved that the density given has negative values at some interval of time if C ; 
If C2 7^ then the Solution III can be written as follows: for % — 1, . . . , n 



Ri{r- V \CI V {T)+K v [ 



T 



-2/[h(D-l)] 



X 



Pi 



21] T 



h 



' CI v+1 (t)-K v+1 (t) 
CI v {t) + K„{t) 

h (D- ly -1-v 
+ 

T 



+ 1 



exp 
1 



are 



-riot 



D- 1 

CI„ +1 (t)-K„ +1 (t) 
CI„{t)+K v {t) 



£ 



CI v+1 (t)-K„ +1 (t) \ 2 
CI u (t)+K v (t) 



1 



(5.29) 
0. 

(5.30) 
(5.31) 

(5.32) 
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2rf D - 2 ™ 



h 2 K 2 D 



j![R- 2Ni T^V \CI v {r)+K v {r)\- h 
1 i=i 



x 



'CI v+1 (t)-K v+1 (t) 
. CUt)+KJt) 



(5.33) 



h 2 n 2 D-2 V \ 1 

1=1 

x f>[C7, +1 (r) - ^ +1 (r)] 2 + [CI„(t) + ^(r)] 2 



(5.34) 

where the independent constants R\, . . . , Rn,C = C\jGi are used. The variable r was 
determined in Eq. (|4.21| ). Also we introduced Kasner parameters e 1 by 



a % 1 
£ l = sgn[/i]— + — 



(5.35) 



By using the properties given in Eq. (4.29) and the asymptotical behavior of the modified 
Bessel functions |TJ, one may prove that the density given in Eq. (5.33) has no negative values 
during the evolution only for : 

1. partial solution with C = on the interval (— oo, +oo) of the harmonic time; 

2. solution with C < on the interval (£o( r o)? +°°)> where r is the root of the equation 

CI v (t o ) + K v (t q ) = 0. (5.36) 



We note that the asymptotical behavior of the solution with C < as t — > +oo is given by 
the partial solution with C — 0. 

Let us now investigate the models having a singularity at the beginning of the evolution. 
Such solutions arise when h < 0. In this case for the solutions with C < one obtain 
the following correspondences (— r) G (— r ,0) £ G (£ (t o ), +oo) £ s G (0, +oo). The 
singularity at £ s = is of Friedman type. The shear viscosity leads to the isotropization at the 
final stage of evolution. We obtain from Eq. (5.30-5.33) exp[x l ] ~ t]^ D ^ and p* = p ~ tj 2 
as t s — > +oo, i.e. the model describes the isotropic Friedman-like expansion corresponding to 
Zeldovich matter (h = 0). The anisotropic behavior is possible on some interval of time after 
the Friedman-like behavior. One can prove that if the constant \C\ ^ is small enough then 
on some interval the sign of the Hubble parameter x l coincides with the sign of the Kasner 
parameter e l (see Eq. (5.31) for h < 0). However, such regime is possible only on limited 
interval of time, because the final stage of the evolution exhibits the isotropic expansion due 
to the shear viscosity. 
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